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Abstract: In this work we calculate the Sunyaev-Zel'dovich (SZ) effect due to the e^e~ 
from dark matter (DM) annihilation in galaxy clusters. Two candidates of DM particle, 
(1) the weakly-interacting massive particle (WIMP) and (2) the light dark matter (LDM) 
are investigated. For each case, we also consider several DM profiles with and without 
central cusp. We generally find smaller signals than previously reported. Moreover, the 
diffusion of electrons and positrons in the galaxy clusters, which was generally thought to be 
negligible, is considered and found to have significant effect on the central electron/positron 
distribution for DM profile with large spatial gradient. We find that the SZ effect from 
WIMP is almost always non-observable, even for the highly cuspy DM profile, and using the 
next generation SZ interferometer such as ALMA. Although the signal of the LDM is much 
larger than that of the WIMP, the final SZ effect is still very small due to the smoothing 
effect of diffusion. Only for the configuration with large central cusp and extremely small 
diffusion effect, the LDM induced SZ effect might have a bit chance of being detected. 
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1. Introduction 

The Dark matter (DM) problem is one of the most important issues in modern physics 
and cosmology. After about eighty years since the first discovery of DM in the Coma 
cluster by Zwicky ||l|, the evidences of DM are overwhelming nowadays. However, most 
of the evidences come from the gravitational effects by astronomical observations, such 
as the rotation curve of spiral galaxies dynamics of galaxy clusters @, gravitational 
lensing effect @, large scale structure of the universe and the anisotropy of cosmic 
microwave background ^. The studies on primordial nucleosynthesis IC] and 



structure formation [^] show that most of the DM is non-baryonic and cold. The nature 
of DM particle is still unknown and remains as one of the biggest puzzles in physics and 
astronomy today. 

Many candidates of DM have been proposed in literature (for reviews, see e.g. |jl2|, [T^]). 
Among the "zoo" of the DM candidates, the weakly- interacting massive particles (WIMP) 
are most favored since they appear natually in many of the new physics models at the 
electroweak scale and can give the correct relic density of DM. The masses of WIMP are 
typically in the range of a few GeV to several TeV, and for weak scale interaction its relic 



density agrees roughly with observation |l^]. The most popular example is neutralino, 



which is the lightest supersymmetric particle in the minimal supersymmetric extension 
of the standard model (MSSM). Typical mass of the neutralino can not be significantly 
lighter than a few GeV However, scenarios with light particles with mass from MeV 
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to GeV are also able to satisfy the constraints from relic density and other astrophysical 
observations 16 1. Furthermore, the light DM (LDM^) annihilation was proposed to 



explain the 511 keV line emission from the bulge around the Galactic center [17|. Further 
studies from astrophysical constraints limited the parameter space of LDM in a narrow 
range |l|, |l|, |20, P, |2|, |3|. 



To identify the nature of DM particles, it is necessary to "see" them in particle physics 
experiments beyond the gravitational measurement. There are usually three types of ex- 
periments suggested to capture the DM particles: the collider-based searches to observe 
the missing energy in particle collisions, the direct searches to capture the scattering sig- 
nals between DM particle and detector nucleus, and the indirect searches to measure the 
annihilation products in cosmic rays like 7-rays, anti-particles and neutrinos etc. (for a 
review see Ref. p^ ). 

It has been suggested that the inverse Compton (IC) scattering between electrons/positrons 
induced by DM annihilation and the CMB photon, which causes a non-thermal Sunyaev- 
Zel'dovich (SZ) effect on CMB, could be an alternative way for DM indirect searches 



1 24, 25, |2g, |2^, |2^. The DM induced SZ effect shows a specific spectrum which differs from 



the usual thermal SZ effect, and can be isolated from other SZ effects with observations of 



arcmin angular resolution and //K sensitivity [28|. In this and a companion paper (paper 
II) we revisit the detectability of the DM induced SZ effect in two typical kinds of objects: 
the galaxy clusters (this paper) and dwarf galaxies (paper II). Both the WIMP and light 
candidate of DM particles are considered. 

The outline of this paper is as follows. In Sec. 2 we introduce the thermal gas distribu- 
tion and DM configuration we adopted, and give the cluster sample we used in this work. 
The production of electrons /positrons from DM annihilation is presented in Sec. 3, and the 
propagation of electrons/positrons in cluster is described in Sec. 4. Sec. 5 gives the results of 
the SZ effect calculations. Finally we draw the conclusions and some discussions in Sec. 6. 
Throughout the paper, we assume a flat ACDM model, with the cosmological parameters 
of the Wilkinson Microwave Anisotropy Probe (WMAP) five year best fitted values, i.e. 
Qm = 0.28, = 0.72 and Hubble constant /i = 0.7 [§. 

2. Gas and DM distributions in cluster 

The galaxy cluster consists of fully ionized gas with temperature 1 ~ 10 keV. The thermal 
electron distribution can be described using an isothermal /3-model [29] 

2 - — 

ne(r) = neo + J|) ' ' (2-1) 

where rigo is the central number density of electrons, rc is the core radius and P represents 
the square of ratio of the galaxy-to-gas velocity dispersions in the cluster |29, 3C]. The 



^Here we give separate discussions about LDM and WIMP according to the usual conventions in the 
Uterature. However, it should be noted that the light DM might be a sub-class of WIMP DM with the only 
difference coming from that it does not suffer from the Lee- Weinberg limit, as emphasized in Refs. p^ . 
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isothermal /3-model is well consistent with the X-ray observations and the parameters can 
be precisely determined from the X-ray surface brightness image of the cluster. 

For the DM distribution in the cluster, however, the most precise knowledge comes 
from numerical simulations. N-body simulations show that there may be a nearly universal 
central cusp of DM density profile, though the exact slope near the center is still being 
debated 32, 33], 34]. However, the observations of the rotation curves of galaxies favored 



cored DM distribution [35, |3^, 37 1. For a review of the DM halo properties please refer to 
Ref. 1 35]. In this work we will adopt the following three types of DM density profiles for 
discussion: 



p{r) = —-^ — , , (hereafter B95, Ref.M), (2.2) 

(l+r/r,)[l + (r/r,)2] ^ ^ ^ 

p{r) = , , I ,2 (hereafter NFW, Ref. P), (2.3) 

(r/rs)(l + rjrsY 

Pi 



p{r) = (hereafter M99, Ref.[g9|). (2.4) 



There are some other profiles with different inner slopes also proposed in literatures ]32 
|33| , Most of these results show similar behaviors (~ r~^) at large radius, but show 
discrepancies in the inner region of the halo. Here we employ B95, NFW and M99 pro- 
files as prototypes of non-cuspy, mediately cuspy and strongly cuspy profiles of DM halos 
respectively. 

Physically, however, the density of DM halo should not diverge, so we introduce a 
cutoff radius within which the density is kept at constant Pmax! probably due to balance 
between the annihilating rate and the in- falling rate of DM J^] . Typically 
10^*^ ~ 10^^ Mq kpc-3 UTI, and we fix /Jmax = 10^^ M© kpc"^ ^.j^jg ^^^.j^^ 



The profile parameters ps and are determined using the virial mass Mvir and con- 



centration parameter Cvir of the DM halo |42]. The virial radius of a DM halo is defined 
as 

1/3 

(2.5) 



,(47r/3)Ap,(z 

where A is the overdensity, and Pc{z) is the critical density of the universe at the redshift 
of the cluster. For ACDM universe, A ISvr^ + 82x - 39a;2 with x = 9.m{z) - 1 = 



— is found to be a good approximation 1 42 ] . The concentration parameter 

is defined as 



' vir 



(2.6) 

d(r^p) 

where r_2 refers to the radius at which ' \r=r_2 = 0. The concentration parameter Cvir 
relates rvir and the density profile parameter as 

^B95 _ ^vir ^NFW _ ^vi£ ^M99 _ ^vir 7) 
1.52Cvir Cvir 0.63Cvir 



Therefore if the Cvir — Afvir relation is specified, is determined using Eq. (|2.7|) , and then 
Ps can be derived by mass condition j p{r)dV = M^i^. Generally the Cvir — M^i^ relation 
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Table 1: The /3-model parameters of thermal electron distribution and total virial mass of the 
cluster sample. 



Name 


z 


neo(cm ^) 


6'c(arcsec)-'^ /3 


Kl [KB\ ) 




Ref.2 




noil 


1 Q V 1 n~3 


441 


0.70 






m 1 


n 


Abell 262 


0.016 


7.25 X 10-3 


87 


0.39 


1.45 


2.70 


m 


471 


Coma 


0.023 


3.42 X 10-3 


624 


0.75 


7.80 


12.9 


m 1 


49| 


Abell 2199 


0.030 


9.90 X 10-3 


132 


0.61 


3.13 


7.00 


m 


471 


Abell 496 


0.033 


5.04 X 10-3 


194.4 


0.64 


6.34 


5.20 


m 1 


El 


Abell 2717 


0.049 


9.60 X 10-3 


35.5 


0.48 


1.64 


1.92 


m 


511 


Abell 1795 


0.063 


2.45 X 10-2 


45 


0.57 


6.74 


9.07 


m 1 




Abell 478 


0.088 


2.36 X 10-2 


50.4 


0.62 


8.32 


16.0 


m 


521 


PKS0745-191 


0.103 


5.50 X 10-2 


24.6 


0.57 


7.70 


11.9 




m 


9c connects with the core radius by the ang 


ular radius dA as Qc ^ 


i rc/dA- 






The former reference 


is for the gas 


distribution and temperature, 


while the latter 


one 





is for the mass. The mass is actually adopted from the compilation of Ref. 1451, where 
the correction of different cosmological models are made. The one listed is the original 
reference. 

should be inferred from the numerical simulations p^, ^ ]. In this work we adopt the fitting 



observational Cvir — Afvir relation 1451 



14.5 / M ^ 



where the reference mass = 1.3 X 10^3/j-i 

We compiled a sample of 9 nearby clusters with redshift z < 0.1 to discuss the non- 
thermal SZ effect from DM annihilation in this work. The /3-model parameters and virial 
mass are given in Table 1. Since the basic properties of the SZ effects are similar for 
different clusters, we will use the Coma cluster as the prime example to introduce the 
calculation process. The results from all of the clusters are presented last. 

3. Electron/positron production from DM annihilation 

The production, propagation and scattering with photons are the same for the electrons 
and the positrons in the current work, so for simplicity we will only discuss electrons in 
the following. The final resulting SZ effect are multiplied by a factor 2 to account for the 
contributions from positrons. The electron source function from DM annihilation can be 
written as 

where {av) is the velocity-weighted annihilation cross section, is the mass of DM 
particle, dA^/d£^ is the electron yield spectrum per annihilation, and p(r) is the density of 
DM. In this work we consider two types of DM particles: WIMP and LDM. 
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3.1 WIMP 



Neutralino is taken as an example of WIMP DM. The direct channel to e^e~ is suppressed 
for neutralino, so electrons are in most cases produced in the cascades of the annihilation 
final-state particles such as heavy leptons, quarks and gauge bosons [13|. The spectra 
of electrons can be different from each other for different annihilation modes. We use the 



package DarkSUSY |54| to calculate the cross section and final-state spectra of electrons. In 
Figure || we plot the electron yield spectra for three typical modes of neutralino annihilation: 
W+W-, hi and r+r". 




E(GeV) 



Figure 1: Electron yield spectra AN/dE for W^W , bb and t+t annihilation modes for neu- 
tralino with mass mwiMP = 100 GeV. 

The annihilation cross section is not well constrained yet. Ref. [^ ] proposed a very 
conservative theoretical constraint {av) < 10~^'^( ^™^y )~^ cm^ s~^. Also there are many 
constraints from the indirect detection experiments |5^, however, all of these re- 

sults depend strongly on the model, including both the DM models and the astrophysical 
background estimations. If the WIMP DM is thermally produced in the early universe. 



the relic density we observe today can provide a fiducial value of the cross section [12| 



3 X 10^27 pj^S g-i 
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with the DM density ^x^"^ = 0.1143 ±0.0034 from the recent combined analysis of WMAP 
five year data together with Type la supernova and baryon acoustic oscillation data 
The DM relic density indicates a cross section (o"v)wimp — 3 x 10~^^ cm^ s^^. Scanning 
the parameter space of MSSM, we can indeed find a series of models which have such a 
cross section and satisfy the relic density condition. For clarity we will fix this value of 
cross section for WIMP DM in the following discussion. 

3.2 LDM 

The annihilation of LDM with mass 1 — 100 MeV was proposed as the source of the 

511 keV line emission in the Galactic center as observed by INTEGRAL |17], though more 

recent observation and analysis seems to favor a more conventional origin of these electrons 

|60| , 61, The constraints from COMPTEL and EGRET observations on the 7-rays 

produced by electromagnetic radiative corrections to XX ~^ e'^e~ process require < 20 

MeV ||l^. More strict mass upper bound of several MeV was found in |21, 22]. Constraints 

from CMB anisotropy yields similar bound [^]. Other bounds include the value of 5 — 2 

|2C], big bang nucleosynthesis (BBN) ||l^. For such low mass DM particles, the allowed 

annihilations can only be to e"*'e~, 7-rays and neutrinos. Following ||T^ we assume that 

the only annihilation channel of LDM is xX ~^ e~^e~. This assumption means that the 

annihilation of the LDM is "invisible" to most other observations ||l^. The final state 

electrons are monochrome with energy spectra dN/dE = 5{E — rriy.)- 

The cross section derived from the flux of the 511 keV 7-ray emission is (o"f)LDM ■ 
2 



IMeV 



10 



-29 



10 cm'^ s 



1 17], which is consistent with the observational relic 



density of DM ]17]. This LDM cross section must be regarded as an upper limit now, as the 
more detailed analysis show that a large fraction of the Galactic center positrons observed 



by the INTEGRAL is produced by astrophysical sources ]^, p2 |. 



4. Propagation of electrons in cluster 

The electrons propagate diffusively in the cluster and experience energy loss processes 
due to the IC scattering, synchrotron radiation. Coulomb collisions and bremsstrahlung 
emission. The propagation equation for electrons can be written as 



V • 



_d_ 

dE 



+ QeiE,r)=0, 



(4.1) 



where due/dE is the equilibrium electron density distribution, D{E, r) and b{E, r) = dE/dt 
are the diffusion coefficient and energy loss rate respectively. For simplicity one can assume 
D{E,r) and b{E,r) to be independent of the spatial location in the cluster. The diffusion 
coefficient can be adopted as a power law with respect to electron energy [^] 



D{E) = Do 



dB \ 
Ikpc J 



2/3 



i^^G E y/^ 

B IGevJ 



(4.2) 



where Dq is a constant, B is the average magnetic field, and ds is the minimum scale of 
uniformity of the magnetic field. For Coma, these parameters are estimated as Dq ~ 3x 10^^ 
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cm^ s ^ , B 1/iG and d-Q ~ 20 kpc [^]. The energy loss rate b{E) is [64] 

m 



10-i6GeV s~^ 



bic{E) + bsjn{E) + 6coui(-E') + "brem 

(E) 



= 0.25 X 
+ 6.13 X 



IGeV 
1 n 



+ 0.0254 X 



(5 lcm~3 
n E \ 



B (3E Y 

1/iG IGeV ) 

( lcm~^ 
1 + 0.013 In 7 



n 



(4.3) 



where /3 and 7 are the velocity and Lorentz factor of electrons respectively, and n is the 
number density of the thermal electrons in the cluster as given in Eg. (2.1). 

4.1 The diffusionless approximation 

In Ref. IMl , it was argued that in galaxy clusters the spatial diffusion is negligible compared 



(4.4) 



with energy loss, then Eq.(4.1) has very simple solution 



drif. 



d^ h{E) Je 



dE'Qe{E',r), 



where Qe{E,r) is the source function given in Eq.(3.1). 




Figure 2: Left: normalized equilibrium momentum (p — 7/?) spectra of electrons generated from 
DM annihilation and thermal electrons with kTe ~ 7.8 keV in Coma cluster. We set ttildm — 5 
MeV for LDM and mwiMP = 100 GeV for WIMP. The three typical channels W+W, bb and 
r+r" for neutralino WIMP annihilation are shown respectively. The central value of the thermal 
electron density is adopted to calculate the energy loss. Right: electron density distributions in 
Coma cluster as functions of radius. Solid line represents the thermal electrons, and the other three 
lines show the DM induced electrons for B95, NFW and M99 profiles respectively. See the text for 
details. 

In the left panel of Figure § we show the normalized equilibrium spectrum f{p) oc ^ 
of energetic electrons from DM annihilation as a function of dimensionless momentum 
p = ji3. The masses are adopted as 100 GeV for WIMP and 5 MeV for LDM. For 
comparison we also plot a thermal electron population with temperature kT^ = 7.8 keV, 
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which is the temperature of thermal electrons in the central region of Coma cluster |48|. 
We see in this figure that at high energies (£■ > 1 GeV, p > 10^), the spectra are very 
soft due to the severe energy loss through the IC scattering. For energy < 100 MeV, the 
Coulomb loss dominates the energy loss processes pi. From Eq.(|43D we know that ftcoul is 



approximately a constant for relativistic electrons, and cx for non-relativistic electrons. 
Therefore we see in Figure ^ that f{p) is flat for 1 < p < 100 and decreases for smaller p. 

In the right panel of Figure |2| we plot the spatial distributions of electrons for both 
the thermal population and DM annihilation induced population in the Coma cluster. 
For WIMP DM the parameters are mwiMP = 100 GeV, (crz;)wiMP = 3 x 10^^^ cm^ 
s~^, and the annihilation channel is W~^W~. For LDM we adopt rra-wiMP = 5 MeV 
and (crt')LDM = 2.5 x 10^^^ cm'^ s^^. It is shown that for thermal electrons the central 
distribution is a flat core, while for NFW and M99 DM scenarios it is strongly cuspied 
in the center. We also note here that at large radius (> 100 kpc), the electron density 
decrease less rapidly than expected from the decrease of DM density(~ r~^) because the 
Coulomb energy loss rate also decreases, thanks to the lower thermal electron density at 
large radii. 

However, even for r < 1 kpc (corresponding to an angular scale < 2"), the density 
of non-thermal electrons produced by DM annihilation is still several orders of magnitude 
lower than the thermal one, therefore it would be very difficult or impossible to detect the 
DM annihilation-induced SZ effect (oc Ue, see Sec. 5). Since the number density of electrons 
is proportional to 1/m'^, the LDM model could generate more electrons, and might be able 
to produce some observational signals of SZ effect. 

4.2 The effect of difTusion 

It is known that the overall diffusion time of electrons in galaxy cluster is much longer than 
the energy loss time ||2^. However, if we focus on the local region such as the central part 
of the cluster, the diffusion term could still be comparable with the energy loss term and 
may significantly modify the equilibrium electron spectrum and spatial distribution. We 
now consider this effect. 

In this case, using the Green's function method, the solution of the electron density 
can be written as 



dUf. 1 



d^ b{E) 

where 



poo 

/ dE'Gir,Av)Qe{E',r), (4.5) 
Je 



1 r' 
G{r,Av) = -j== (-1)" / dr'- x 

{r'-rnY\ ( (r' + r„) 



exp — — exp 



4At; / V 4A?; 



^2 



(4.6) 



r„ = (— l)"r -|- 2nrh is the location of the nth "charge" image, is the radius of the 
diffusion halo, and Av{E,E') = f^'deD{e)/b{e). 

In the left panel of Figure ^, the momentum spectra of the electrons with the diffusion 
effect included are shown as lines. For comparison, we also plot the results without diffusion 



-8- 



10 10"' 10° 10' 10^ 10' 10* 10"' 10° io' lo'^ 10° 

P r(kpc) 



Figure 3: Left: the normalized electron spectra f{p) for r — 0.1, 10 and 10"^ kpc respectively, 
compared with the results without diffusion (shown by points with the same colors as the lines). 
Right: the electron density profiles for three energies E = 0.1, 1 and 50 GeV, compared with the 
results without diffusion (points). In this calculation we adopt WIMP DM with W^W~ channel 
and M99 density profile. See the text for details. 

as points on the same figure. In this calculation we adopt a WIMP DM with mass = 
100 GeV, assuming the W^W~ channel, and the M99 density profile. We see that the 
diffusion indeed can lead to distortion of electron spectrum, especially when r is small. For 
r > 10 kpc, where the spatial gradient of the DM distribution is small, the effect becomes 
negligible. The spatial distributions of the electrons for several energies are shown in the 
right panel of Figure |3|. It is clear that at small radii diffusion leads to a smoothing of the 
spatial profile of electrons, while at large radii the results approach the diffusionless solution. 
We also note that the differences between the cases with and without diffusion are energy- 
dependent. For higher energy, the energy loss is more important, and the diffusion effect is 
weaker, so the differences begin to appear at smaller radii for higher energy. An important 
consequence is that the strong central cusp of electron distribution is smeared out, which 
significantly affects the SZ effect for high angular resolution observation. Similarly, for the 
NFW profile there is also a smoothing effect from diffusion. Only in the case of the cored 
B95 profile the diffusion does not make a difference. 



where x = hv/kT^ is the dimensionless frequency of CMB photon, Tq = 2.725 K is the 
undistorted CMB temperature, g{x) is the spectral distortion function, and y{0) is the 
Comptonization parameter which is proportional to the number density of electrons, for 
angle separation 9 from the center. The spectral distortion function can be expressed as 



5. SZ effect 



For the calculation of SZ effect, we follow the method presented in Ref. The temper- 
ature variation of CMB after traveling through a population of electrons is 





rrif-C 



.2 r 




g{x) 



io{xe ^)Pi(s)ds — io(x) 
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where io{x) = x^/{e^ — 1) is the Plankian distribution of CMB photons, s = ln(z^'/z^) is the 
frequency shift of one photon after one scattering with electrons, and Pi{s) is the frequency 
shift probabihty distribution, 

Pi{s) = I f{p)Ps{s,p)dp, (5.3) 

with p and f{p) are the dimensionless momentum and momentum spectrum respectively, 
and Ps{s,p) is the probability of a photon has an energy shift s when colliding with an 
electron with momentum p j6^, The (kTe) in Eq.( |5.2D is the average effective electron 



temperature defined as |65] 



(fcTe) = = / -f{p)pPm,c'dp, (5.4) 

J Hedl J S 

with Pe the pressure of electrons, and dl the line-of-sight integral. Finally the Comptoniza- 
tion parameter y(6) in Eq.(|5.1|) is given by 



y(e) = = ^.aT [ nedl, (5.5) 

where aj- is the Thomson cross section. For the case with diffusion effect, the spatial and 
energy distributions of electrons are coupled together, g{x)y{9) is integrated together in 
the line-of-sight integral. 

Since the DM induced electrons will concentrate near the center of the cluster for cuspy 
density profiles, while from the observational point of view we only have limited resolution 
angle, it is necessary to smooth the results within the resolution angle of the detector array. 
In Figure Q we show the average Comptonization parameter, ?/sm(^sm) = jSsm g^g — > ^ 
function of the beam size angle ^sm for the WIMP DM. It is shown that for the diffusionless 
case ysm varies approximately with 9^^, and for B95, NFW and M99 profiles 
respectively. For the case with diffusion the central cusp of electrons from NFW and M99 
profiles are smoothed out and the angle dependence of ysm becomes much weaker. 

The expected temperature variation of CMB due to the e~^e~ from DM annihilation 
in Coma cluster is shown in Figure |5[ The left panel is for WIMP with mass ttt-wimp = 
100 GeV and annihilation channel W~^W~; the right panel is for LDM with mLDM = 5 
MeV. The beam width is assumed to be 1 arcsec around the center of the cluster, which 
corresponds to a radial distance ~ 0.5 kpc from the center. In each panel, we plot the 
thermal SZ effect as a solid red curve (on the left panel, due to the scale of the plot, it 
appears as almost a verticle line). The DM-induced SZ effect for the M99 profile without 
inclusion of diffusion is shown as open circles, the DM-induced SZ effect for the M99 profile 
with diffusion as blue dashed curves. The effect for the NFW and B95 profiles are also 
plotted, though as they are much smaller compared with the more cuspy M99 profile, the 
two curves almost coincides with the x-axis and are hardly visible. 

We can see that for typical WIMPs (left panel of Fig. |5|) the DM-induced SZ effect is 
extremely small. At 217 GHz (x ~ 3.83) where the thermal SZ effect is zero, even without 
including the effect of diffusion, the DM-induced SZ effect for the M99 profile is only 
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Figure 4: Smoothed Comptonization parameter for WIMP as a function of the smooth angle 6*8111 
for B95, NFW and M99 density profiles respectively. The thick lines show the results without 
diffusion effect, and the thin lines are the cases with diffusion effect. 
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Figure 5: CMB temperature variation due to the e~^e~ from DM annihilation in Coma cluster. 
Left panel is for WIMP with mass mwiMP = 100 GeV (VF+Vl^" channel), and right panel is for 
LDM with rnLDM = 5 MeV. See text for detail. 



-3 X 10 '^fiK. These results are much smaller than those given in pi, E6|^. While we have 



^Note that in a recent study similar conclusion is also derived 
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adoptted different parameter values in the the models presented above, we have also checked 
the cases given in i.e., N04 density profile = 40 GeV, av = 4.7 x 10~^^ cm^ s~^ 

and bb annihilation channel. We find that the temperature variation to be 
for frequency 30GHz, which is about three orders of magnitude lower than the result of 
~ 40^K given in Ref.@. 

For the case of LDM (see right panel of Fig. ^), the DM induced SZ effect is larger. 
In the diffusionless approximation, at frequency 217 GHz, for example, the temperature 
deviations for M99 profile is — 16/uK, which is comparable with the sensitivity of the next 
generation SZ interferometer such as the Atacama Large Millimeter Array (ALMA)^. 

However, with the effect of diffusion which lowers the density of electron-positron pairs 
near the center of the cluster, the DM-induced SZ effect is even smaller. For WIMP, even 
for the M99 profile, it is only -2 x 10~VK (about IQ-^/xK for NFW profile, lO^^^uK for 
B95 profile). For LDM, although it is larger than the WIMP case, the DM-induced SZ 
effect is only — 0.75/iK for the M99 profile, making its detection extremely difficult if not 
impossible. 

Finally we plot in Figure ^ the calculated SZ effects for the 9 clusters given in Table 
1. The DM profile is adopted as M99 and the results are smoothed within 1 arcsec around 
the center of the halo. Note that the non-thermal SZ effect from DM annihilation is 
more remarkable for nearby clusters, because for far away clusters the same beam size will 
correspond to a larger region around the center, and the average effect becomes smaller. 
We can see that if the diffusion effects are taken into account, there is almost no chance 
to detect the DM induced SZ effects due to the very weak signals and strong thermal 
backgrounds. 



6. Discussion 

In this work we calculate the SZ effect induced by e^e~ from DM annihilation in galaxy 
clusters. Two types of DM particles, WIMP and LDM, are considered. The annihilation 
cross sections we adopted satisfy the constraint from the relic density of DM. Neutralino in 
the framework of supersymmetry is taken as an example of WIMP DM, and three typical 
annihilation channels, W^W~ , bb and t~^t~ are employed. For LDM we assume the only 
annihilation channel is e~^e~ , and its annihilation cross cross to produce all of the positrons 
in the Galactic center, which should be regarded as very conservative upper limits. The 
density profile of DM halo is a crucial factor, and we consider the Burkert, NFW and 
Moore profiles to represent the case of non-cuspy, mediately cuspy and strongly cuspy 
profiles respectively. 

We find much smaller (two orders of magnitude or more) DM-induced SZ effect than 
previous claims for WIMPs. Furthermore, we consider the spatial diffusion of electrons, 
which was negelected in previous works. This effect significantly reduced the density of 
energetic electrons in the center of halo, especially for the DM profile with strong cusp. 
Due to this effect, the DM induced SZ effect is even weaker. For WIMPs, the DM induced 

^http:/ /www. eso.org/sci/facilities/alma/iridex. html 
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Figure 6: Temperature variations of CMB due to tlie e+e from DM annihilation for the 9 clusters 
listed in Table 1. The masses of WIMP and LDM are 100 GeV and 5 MeV respectively, and the 
DM profile is M99. 



SZ effect is several orders of magnitude smaller than the thermal one. Although the LDM 
could produce stronger signal than the WIMPs, it is still very small. Given the small size 
of the DM-induced SZ signal, not to mention the practical issues of seperating the signal 
from possible confusions such as astrophysical foreground, kinetic SZ effect, and fluctua- 
tion of primordial CMB temperature, we conclude that it would be extremely difficult if 
not entirely impossible to detect the DM induced SZ effect with the next generation SZ 
interferometer such as ALMA. 

There are large uncertainties about the particle physics properties of DM such as the 
mass and annihilation cross section. For WIMP, we fix the cross section to be (o"t;)wiMP — 
3 X 10~^^ cm^ s~^ taking into account the constraint from the relic density, so the SZ effect 
will be approximately proportional to l/m^jj^p. For lower mass WIMPs, the SZ effect can 
be stronger. Furthermore, larger cross section is also possible in some scenarios such as the 



non-thermal production of DM or the "Sommerfeld enhancement" |6£, 70, 71, 73|, but 



these scenarios is constrained by, e.g., the 7-rays [74, 75]. For LDM, since the annihilation 



cross section derived from the 511 keV observations at the Galactic center scales as w-ldm, 
the change of electron density through varying mlj^M compensated by a rescaled cross 
section. However, different mass will result in different high energy cut off for the equi- 
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librium electron spectrum (Eq. (|4.4| )) and the eventual SZ effect would be affected^. More 
quantitatively we find that the temperature distortion is nearly proportional to rriLDM- 
Since the mass of LDM is constrained in a narrow range by other methods |18, |2^, 21, 23|, 
we think the basic conclusion of this work is still valid. 



The DM substructures inside the halo might also "boost" the annihilation signal |76 



77| , |78| . However, the substructures near the center of the halo would be destroyed by the 



tidal force, and the "boost" can only take effect at large radii [^. Furthermore, for the 
cuspy profiles the central density is high enough to dominate the contribution over that 
from substructures, hence the "boost" effect is relatively weak. Therefore if we investigate 
the SZ effect at the center of the cluster, the effect of substructures can be reasonably 
neglected. 

Finally we note that the recent observations of the electrons/positrons by ATIC and 
PAMELA experiments show apparent excesses of energetic electron/positrons compared 
with the conventional cosmic ray model predictions []80| , |8l|| . If these excesses are ascribed 
to DM annihilation, the mass of DM particle > 700 GeV and a boost factor of order ~ 10^ 
are needed, and the annihilation modes should be lepton dominated (e.g., Ref. ^2[). We 
also calculate the SZ effect of this DM configuration, and find that the results are almost 
the same as the 100 GeV WIMP used in this work. This is mainly because the boost factor 
is canceled by the 1/m^ term. That is to say, the SZ effect from DM annihilation is still 
very difficult to be detected even for the enhanced electron/positron density as required 
by ATIC and PAMELA. 
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